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A PIECEWISE LINEAR FEM FOR AN OPTIMAL CONTROL 
PROBLEM OF FRACTIONAL OPERATORS: ERROR ANALYSIS 

ON CURVED DOMAINS 

ENRIQUE OTAROLA 


Abstract. We propose and analyze a new discretization technique for a linear- 
quadratic optimal control problem involving the fractional powers of a symme¬ 
tric and uniformly elliptic second oder operator; control constraints are con¬ 
sidered. Since these fractional operators can be realized as the Dirichlet-to- 
Neumann map for a nonuniformly elliptic equation, we recast our problem 
as a nonuniformly elliptic optimal control problem. The rapid decay of the 
solution to this problem suggests a truncation that is suitable for numerical 
approximation. We propose a fully discrete scheme that is based on piecewise 
linear functions on quasi-uniform meshes to approximate the optimal control 
and first-degree tensor product functions on anisotropic meshes for the opti¬ 
mal state variable. We provide an a priori error analysis that relies on derived 
Holder and Sobolev regularity estimates for the optimal variables and error 
estimates for an scheme that approximates fractional diffusion on curved do¬ 
mains; the latter being an extension of previous available results. The analysis 
is valid in any dimension. We conclude by presenting some numerical experi¬ 
ments that validate the derived error estimates. 


1. Introduction 

In this work we shall be interested in the analysis of a new solution technique for 
a linear-quadratic optimal control problem involving the fractional powers of an 
uniformly elliptic second oder operator; control constraints are considered. Before 
describing such a discretization technique, we follow [Sj and illustrate the PDE- 
constrained optimization problem we are interested in. Let f 1 be an open and 
bounded domain in (n > 1) with Lipschitz boundary dQ. Given s £ (0,1) and 
a desired state lq : fl —> R, we define the cost functional 

(1-1) ^(U: z ) = 2 II u — Ud lli 2 (n) + 2 H z lli 2 (fip 

where i? > 0 denotes the so-called regularization parameter. We shall be concerned 
with the following PDE-constrained optimization problem: Find 

(1.2) min J(u,z), 
subject to the fractional state equation 

(1.3) £ s u = z in fl, u = 0 on dil, 
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and the control constraints 


(1.4) a < z(x') < b a.e. x' £ tt. 

The bounds a, b £ R. satisfy the property a < b. The operator C s , with s £ (0,1), is 
a fractional power of the second order, symmetric and uniformly elliptic operator 

(1.5) Cw = —<l\v x ’{A{x')'S7 x 'w) 


supplemented with homogeneous Dirichlet boundary conditions. The matrix of 
coefficients A{x') is bounded and measurable in Q and satisfies a uniform ellipticity 
condition. For convenience, we will refer to the optimal control problem (1.2)- (1.4) 
as the fractional optimal control problem. 

Recently, the analysis of regularity properties of the solution to fractional PDE 
has received a tremendous attention: fractional diffusion has been one of the most 
studied topics in the past decade Qu m m ini nsi- Concerning applications, 
capturing the essential behavior of fractional diffusion is fundamental in science 
and engineering since it allows for understanding of applications where anomalous 
diffusion is considered [Tj, complex phenomena in mechanics [Ti, turbulence |16j . 
nonlocal electrostatics j33]J, image processing [26] and finance [35]. In many of these 
applications, control problems arise naturally M- 

Regarding the numerical analysis literature, the approaches advocated for frac¬ 
tional diffusion are mainly divided in two categories. The first one is based on an 
approximated spectral representation that requires the solution of a large eigen¬ 
value problem [311 |32l 09]. The second approach is based on the singular integral 
definition of fractional diffusion: the obtained matrices are dense and the presence 
of a non-integrable kernel requires special attention [3 03 S3]- As opposed to these 
approaches and inspired in the breakthrough by L. Caffarelli and L. Silvestre m, 
the authors of m have designed and analyzed a new solution technique for frac¬ 
tional diffusion via an extension problem. Although the analysis of the proposed 
method is intricate, its implementation is done using standard components of finite 
element analysis [40] . 

In contrast, the analysis of optimal control problems involving fractional diffusion 
and nonlocal operators is still in its infancy; see ISJEUHIM for some preliminary 
results. In [19] . a PDE-constrained optimization problem involving a nonlocal diffu¬ 
sion equation is considered and analyzed. For such a problem, the authors propose 
a numerical scheme to approximate its solution and derive convergence results. In 
0, the authors have used the Caffarelli-Silvestre extension \TT\, and its develop¬ 
ments to both bounded domains [13 IHj and a general class of elliptic operators 
[131 146] , to design and analyze solution techniques for the fractional optimal control 


problem (1.2) (1.4). The analysis hinges on the fact that the fractional powers of 
C can be realized as an operator that maps a Dirichlet boundary condition to a 
Neumann condition via an extension problem on C = fl x (0, oo). This extension is 
the following local boundary value problem (see [13 [EE] Q. 03 for details): 


(1.6) CW — — dyytf/ = 0 in C , ^ = 0 on OlC, dfW = d s z on f l x {0}, 

where OlC = dfl x [0,oo) is the lateral boundary of C, a = 1 — 2s £ (—1,1), 
d s = 2°T(1 — s)/r(s) and the conormal exterior derivative of % at fl x {0} is 

(1.7) = - lim y a %. 

v- >o+ 
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The limit in (1.71 is understood in the distributional sense mmm • We will call 
y the extended variable and call the dimension n+1 in R" +1 the extended dimension 


of problem (1.6). As noted in ns m El E5], the fractional powers of the operator 
C and the Dirichlet-to-Neumann operator of problem (1.6) are related by 


d*C s u = 3“^ in 1 1. 


Note that the differential operator in (1.6) is —div (j/“AV^), where A (x',y) = 
diag{A(a/), 1} for all (x r , y) in C. Thus, we rewrite (1.6) as 


(1.8) -div(y Q AV^) =0 in C, ?/= 0 on d L C, 3“^ = d s z on SI x {0}. 


Invoking such a localization result, the authors of [5] overcome the nonlocality 
of C s by considering an equivalent formulation for ( |1.2| )-( 1.41: min J(^| y=0 ,z) 
subject to the linear equation (11.8b and (1.4); the extended optimal control problem. 


On the basis of this formulation, the following simple strategy to solve (1.2) (1.4) 
is proposed: given s £ (0,1) and a desired state : il —► R, solve the extended 
optimal control problem, thus obtaining z(x') and &(x',y). Setting 0 : i' € !1 £ 
u(x') = ^(x',0) £ R, the optimal pair (u,z) solving ( |1.2| )-( 1.4) is obtained. Two 


numerical techniques are analyzed in [5] : one that is semidiscrete where the optimal 
control is not discretized, and the other one that is fully discrete and discretizes 
the optimal control using piecewise constant functions. A priori error estimates 
are derived; the ones for the fully-discrete scheme being quasi-optimal in terms of 
approximation. 

In this work, we continue with our recent research program by proposing and 
analyzing a new solution technique for the optimal control problem (fW2T) (p^4l). 
This technique corresponds to a fully discrete scheme, which, in contrast to 0, 
discretizes the optimal control with piecewise linear functions on quasi-uniform 
meshes. Following £40], the optimal state is discretized with first-degree tensor 
product finite elements on anisotropic meshes. At this point, it is important to 
comment that although the essential localization results have been already intro¬ 
duced and analyzed in [5], the error analysis of our proposed scheme comes with 
its own set of difficulties. Overcoming them has required us to provide several new, 
nontrivial, results. Let us briefly detail some of them: 


1. We provide regularity results, in Holder and Sobolev spaces, for the continuous 
optimal control variables that require some suitable assumptions on the smooth¬ 
ness of the domain and the matrix of coefficients A. 

2. In order to invoke the regularity results detailed in the previous point, we have 
seen the need of developing, inspired in j40| . an a priori error analysis for frac¬ 
tional diffusion on curved domains f l of class C 2 . The analysis is able to deal 
with both the natural anisotropy of the mesh in the extended variable y and the 
nonuniform coefficient y a with a £ (—1,1). 

3. We provide an a priori error analysis for our scheme that requires a subtle 
interplay between Sobolev and Holder regularity of the optimal control variables, 
as well as as growth conditions, and relies on the assumption that the boundary 
of the active sets consists of a finite number of rectifiable curves 13Z1I3H]. 

It is thus instructive to compare the error estimate derived in this work with the 
one obtained in [5J Corollary 5.17]. The aforementioned error estimates read 


\\z-Z\\ L 2 {n) <\logN\ 2s N-^, and ||z - Z\\ L 2 (n) < | log N\ 2s N~M +e ), 
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respectively. Here z denotes the optimal control, Z denotes the corresponding 
control approximation and N denotes the number of degrees of freedom of the 
underlying mesh. For s £ (1/2,1), the parameter cr corresponds to 1. Meanwhile, 
for s £ (1/4,1/2), we have that a = 2s. Consequently, for s £ (1/4,1) our proposed 
scheme over-perform the one developed in [5]; see Section [7] for an illustration. For 
s £ (0,1/4], we have that cr = 1/2 and, for these values of s, we derive a linear 
order of convergence for the proposed scheme. 

The outline of this paper is as follows. In Section [2] we recall the definition of 
the fractional powers of elliptic operators via spectral theory and introduce some 
suitable functional framework. In Section [3] we define the fractional and extended 
optimal control problems. Section |3.1| presents regularity results for the optimal 
control variables, in both Holder and Sobolev spaces. Section [4] is advocated to the 
truncated optimal control problem. We also state the approximation and regularity 
properties of its solution. In Section [5] we review the a priori error analysis of [JD. 
for the state equation (1.8) and in Section 5.1 we extend these results to a scenario 
where fl is a C 2 curved domain. In Section [b] we propose a fully-discrete scheme 
for the fractional optimal control problem that discretizes the optimal control with 
piecewise linear functions. For s £ (0,1), we provide an a priori error analysis for 
the proposed scheme. The analysis is valid in any dimension. Finally, in Section 
[71 we present numerical experiments that illustrate the a priori error analysis of 
Section © 


2. Notation and preliminaries 

2.1. Notation. Throughout this work Q is an open, bounded and connected do¬ 
main of R n (n > 1) with Lipschitz boundary dtt. We will follow the notation of 
13 HU and define the semi-infinite cylinder with base H and its lateral boundary, 
by C = 12 x (0, oo) and 8lC = dfl x [0, oo), respectively. Given <y > 0, we define 
the truncated cylinder C<y = H x (0,fT) and its lateral boundary di,Cy accordingly. 

Since we will be dealing with objects defined in R" +1 and the extended n + 1 
dimension will play a special role, we denote a vector x £ R n+1 by 

x = (x 1 ,...,x n , x n+1 ) = (x', x n+1 ) = (x\ y), 

with Xi £ R for i = 1,..., n + 1, x' £ R” and y £ R. 

The matrix of coefficients A(x'), that defines the operator C in 
able and bounded in il, and uniformly elliptic. We denote by C s , with s £ (0,1), 
the fractional powers of the operator C. The parameter a belongs to (—1,1) and 
is related to the power s of the fractional operator C s by the formula a = 1 — 2s. 

Finally, the relation a < b indicates that a < Cb, with a constant C that depends 
neither on a or 6 nor the discretization parameters. The value of C might change 
at each occurrence. 

2.2. Fractional powers of a second order elliptic operator. We consider the 
definition based on spectral theory [151 Hil HO] . The operator C 1 : L 2 (Q) —¥ T 2 (H), 
that solves Cw = f in H and w = 0 on <9f2, is compact, symmetric, and positive, 
so its spectrum {A/' 1 }fc g N is discrete, real, and positive and accumulates at zero. 
Moreover, the eigenfunctions 

(2.1) {^fclfceN : £Vk = A fe ( p k in fl, <Pk = 0 on d Cl, k £ N, 


(1.5), is measur- 
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form an orthonormal basis of L 2 (fl). With this spectral decomposition at hand, 
the fractional powers of C can be defined, for w G Cq°(D), by 

OO 

(2.2) C s w = ^2 K w kVk, 

k =1 

where u>k = f n wtpk and s G (0,1). By density C s can be extended to the space 

{ OO OO 

w = ^ Wk<Pk ■ ^2 AfcWfc < oo 
k —1 fe=l 

which, as a consequence of the theory of Hilbert scales [M., coincides with the 
space [Hq (H), L 2 (fl)]i_ s . If s G (1,2], and f l is, for instance, an open bounded and 
convex domain, we have that HF(Q) = H s (fl) n .ffg(fl); see [25] ■ For s G (0,1) we 
denote by the dual space of H s (fl). 


2.3. The nonuniformly extension problem. The localization results of Caf- 
farelli and Silvestre uni HUE] require us to deal with the nonuniformly elliptic 
problem (1.81 and consequently with Lebesgue and Sobolev spaces with the weight 
y a for a G (—1,1). Following (2TJ [48], we define, for an open set D C R n+1 , the 
weighted Lebesgue space 


7 D) = IwG L\ oc (D) : \\w\ 


li 2 (|y| c 


,D) 


= / \y\ a w 2 < oo 


ID 


and the weighted Sobolev space 

H\\y\ a ,D) = {wG L 2 (\y\ a , D) : |Mliyi(|i/|“,D) < oo} , 

where 

(2-4) |M|i? 1 (|l/|“,.D) = (lMli 2 (|j/|“,.D) + l|Vw||l 2 (|3/|“,£>)) 2 • 

Since a = 1 — 2s G (—1,1), we have that \y\ a belongs to the Muckenhoupt 
class A 2 (K" +1 ) EUEl ESS SHI- Consequently, Lf 1 (|i/|“, D) endowed with the 
norm ( |2.4| ) is Hilbert and C°°(D ) fl H l (\y\ a : D) is dense in lL 1 (|r/| a , D) (cf. [48] 
Proposition 2.1.2, Corollary 2.1.6], [34] and [27]. Theorem 1]). For convenience, we 
recall the definition of the Muckenhoupt class A 2 i|39] HE]. 

Definition 2.1 (Muckenhoupt class A 2 ). Let w be a weight, that is cu G Ll oc (R n ) 
and uj > 0 a.e. in R n with n > 1. We say that to G H 2 (R n ) if 

(2.5) C 2 ^ = sup (J- lu'J (J- oj- 1 ^ < oo, 

where the supremum is taken over all balls B in K n . 


We now define the suitable weighted Sobolev space to analyze problem (1.81: 
#i(l/“,C) = {wG H\y a ,C) : w = 0 on d L C} . 

On this space, we have the following weighted Poincare inequality 


(2-6) IMU 2 ( y“,C) l|Vu)|| L 2 ( . yQjC) Vw e Hl(y a ,C); 
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see HOI inequality (2.21)]. Consequently, the seminorm on H\[y a ,C) is equivalent 
to the norm (2.4). For w £ H l (y a ,C), tr^ w denotes its trace onto fl x {0}. We 
recall that, for a = 1 — 2s, EDI Proposition 2.5] yields 

(2.7) tin Al(y a ,C) =H s (ft), || tr n w|| H »(n) < Ctr n |MI#i (l ,a )C) - 


The weak formulation of problem (1.8) reads as follows: Find % £ H\(y a ,C) 
such that 

(2.8) a(°tt, <ff) = (z, trn (>)i-«(n)xi'(0) £ iJ^(y“,C), 

where, for w,(f> £ Hfyy a ,C), the bilinear form a is defined by 

(2.9) a{w, (j>) = f j/ a A(i)Vw ■ V^i, 

d s Jc 

and (■, ■)i-.(n)xiqn) denotes the duality pairing between H s (f2) and H _s (f2) which, 
as a consequence of (2.7), is well defined. We recall that the matrix of coefficients 
A(x',y) = diag-j^a;'), 1} for all ( x',y) £ C. 

The conditions on A and (2.6) imply that a is bounded and coercive in H]^{y a ,C). 
Consequently, the well-posedness of (2.8) follows from the Lax-Milgram Lemma. 
In the manuscript, we will use repeatedly that a(w,w) 1//2 is a norm equivalent to 
(2.4). We present the following estimate for problem Proposition 2.1]: 

(2-10) < ||u|| H «(n) < ||z||h—( n)- 

Define the Dirichlet-to-Neumann operator N : H s (f2) — > H _s (fl) by 
u £ H s (12) i—s> N(u) = £ IT s (f2). 

We thus present the Caffarelli-Silvestre extension result mmm- 

Theorem 2.2 (Caffarelli-Silvestre extension result). If £ H\(y a ,C) and u £ 
HP(f2) solve (2.8) and (1.3) for z £ H - ^!!), respectively, then u = tr^ ^ and 

d s C s u= lV(u) =<9“^ in 

3. The fractional and extended optimal control problems 


For J defined in (1.1), the fractional optimal control problem reads as follows: 
Find min J(u,z) subject to the fractional state equation (1.3) and the control con¬ 
straints & The set of admissible controls Z a( j is defined by 

Z ad = { w £ L 2 (fl) : a < w(x') < b, a.e. x 1 £ f1}, 

where a, b £ R. and satisfy a < b. The desired state u d £ L 2 {fl) and d > 0. 

We define the fractional control to-state operator map S : H -S (fl) 9 z u £ 
H s (f2), where u = u(z) solves ( |1.3[ ). The operator S is linear and, as a consequence 
of (2.10), bounded from H _s (fl) into 10^(12). Given a control z £ H _S (D), we define 
the fractional adjoint state p = p(z) £ ]HF(f2) as p = S(u —u d ). With these elements 
at hand, we recall the following result 0 Theorem 3.4]. 

Theorem 3.1 (existence, uniqueness and optimality conditions). The optimal con¬ 
trol problem ~( |1.4[ ) has a unique optimal solution (u,z) £ H s (f2) x Z od . The 
optimality conditions u = Sz £ H S (S2), p = S(u — u d ) £ H S (I2) ; and 

(3.1) z£Z ad, (i9z+p,z-z) L 2 (n) > 0 Vz £ Z ad 

hold. These conditions are necessary and sufficient. 
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Since the fractional powers of C are non-local operators, the design and analysis 


of solution techniques for problem (1.21 (1.41 is far for being trivial. To circumvent 


such a difficulty, in reference [5] the authors invoke the localization results stated in 
Theorem |2.2| and propose the following equivalent optimal control problem: Find 
min J(tro ^,z) subject to the linear and extended state equation 

(3.2) a(fy,(f>) = (z,tr n i)i)i-. ( n )x i. ( t! ) £ #i(y“, C) 


and the control constraints z £ Z a( j; the bilinear form a is defined as in (2.9). We will 


refer to this problem as the extended optimal control problem. We recall that the 
trace operator tro is defined in Section 2.3| The equivalence of the aforementioned 
optimal control problems follows from [HI Theorem 3.12]: u(z) = tro ^(z). 

We define the extended control-to-state operator G : H _s (f2) 9z£ tro £ 
H s (f2), where = ^(z) solves (|3.2[). Invoking Theorem 2.2 we conclude that the 


actions of G and S coincide and then that G is a linear and continuous mapping 
from into EF(f2). 

We also define the extended adjoint state 2A = ^(z) £ H\{y a ,C), associated to 
the extended state = ^(z), as the unique solution to 

(3.3) a(<p, &) = (tr n ^ - u d , tr n </>) L 2 (n) V</> £ Hl(y a ,C). 

With this definition at hand, we present the following result [5] Theorem 3.11]. 

Theorem 3.2 (existence, uniqueness and optimality system). The extended opti¬ 
mal control problem has a unique optimal solution (^,z). The optimality system 

' & = ^(z) £ H\(y a ,C) solution to (3.2), 


& = &(z) £ Hl(y a ,C) solution to ( |3.3[ ), 

.zeZ ad, (tro & + $z, z — z) L 2 ( q) > 0 Vz £ Z ad , 

holds. These conditions are necessary and sufficient. 

We conclude this section with a representation formula for ^ that is based on 


the eigenpairs {A k,Tk} defined in (2.1). Let u(x') = Yhk u kTk{x') be the solution 
to 11.3|) , then the solution % to problem (13.2|) can be written as follows: 


(3.4) 


OO 


W(x,t) = ) u k ip k (x')if k {y), 


k =1 


OO. 


where, for G N, ^ solves 

(3.5) ipk+ay^ip'k-\ k ip k = 0, ip k ( 0) = 1, tpk(y) -t 0 as y 
If s = 1/2, then ip k (y) = e~'^ v . If s £ (0,1) \ {1/2} and c s = 2 1_s /r(s), then 

f>k{y) = c s (v^y) K s (^/%y), 

where K s is the modified Bessel function of the second kind [3 Chapter 9.6]. We 
refer the reader to m Section 2.4] and [TT] Proposition 2.1] for details. 

3.1. Regularity of the fractional optimal control. Since we will be concerned 
with the approximation of the solution to the extended optimal control problem, it 
is essential to study the regularity properties of (u, p, z). To derive such results, we 
will assume different conditions on the smoothness of the domain 12 and the matrix 


of coefficients A defining C in (1.51; we will be precise regarding these assumptions 


when is needed. We start by recalling the results of Lemma 3.5 in [5]. 
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Lemma 3.3 (regularity of the optimal control). Let z £ Z be the fractional 
optimal control. If LI is convex, A £ C' 0,1 (f2). u d £ H 1-S (f2) and a < 0 < b, then 
z £ 


We recall the so-called projection formula for the optimal control z. If i9 > 0 
and p denotes the fractional optimal adjoint state, then the projection formula 


(3.6) 


z(x') = proj [a b] (-^P(®0 


is equivalent to (3.1), where projj a b ; (v) = minjb,max{a,w}}; see [47l Section 2.8]. 
We comment on the assumption a < 0 < b of Lemma |3.3[ the bootstrap argument 


developed in the proof of Lemma 3.5 in [5] relies on the projection formula (3.6) 
and thus requires the assumption a < 0 < b in order to preserve the boundary 
values of z £ H s (fl). We refer the reader to the proof of j5j Lemma 3.5] for details. 

On the basis of the projection formula (3.6) we derive the following regularity 
properties, on Sobolev spaces, for z and p. 


Lemma 3.4 (regularity of the optimal state and adjoint state). If u d £ H 1_s (fl), 
LI is convex, a < 0 < b and A £ C 0,1 (fi), then u £ H K (fi) for k = minjl + 2s, 2}. 
In addition, if Ud £ Hq(LI), then p £ H K (fl) for k = minjl + 2s, 2}. 


Proof. An application of Lemma 3.3 yields z £ Hq(LI). Then, since Ll is convex, we 
invoke the fact that C s is a pseudodifferential operator of order 2s and Theorem 
3.2.1.2 in [21], to conclude that u, the solution to problem (1.3), belongs to H K (S7), 
where n = minjl + 2s,2}. If u d £ Hq(LI), we have that u — Ud £ .Hq (12) an d then a 
similar argument shows that the optimal adjoint state p £ H K (f2). □ 


In Section [6] we will propose a numerical scheme to approximate the solution to 


(1.2) (1.4) and also provide an priori error analysis. The latter is based on Holder 
regularity results for the optimal control z. In order to obtain such results, we 
invoke the recent boundary regularity estimates for the solution u to problem (1.3), 
derived by Caffarelli and Stinga in ]l3], that assume some suitable smoothness 
properties on z, the domain Ll and the matrix A. Since, in our setting, we have 
that z £ Z ad C L°°(Ll), these results can be adapted to the following regularity 
estimates; see also [T31 Lemma 2.10] and [8], Proposition 5.2]. 


Lemma 3.5 (boundary regularity for u). Let u be the solution to ( |1.3| ) with z £ 
L°°(Ll). If s £ (0, |), Ll is a C 1 domain and A £ C(Ll), then u £ C°’ 2s (Ll), and 

(3- 7 ) [u] C o, 2 S (n) < ||u|| H »(n) + ||z||i,°°(r 2 )- 

On the other hand, if s £ (|, 1), Ll is a C 1 ’ 25-1 domain and A £ C'°’ 2s_1 (f2), then 
u£C 1 ’ 2s ~ 1 (Ll) and 


[ u ]cu 2 »-i(n) ^ IM|H»(n) + ||z|U”(f 2 )- 

When s = \, it is not completely evident how to adapt the techniques developed 
in m to derive regularity results for the solution u under the assumption that 
z £ L°°(Li). For this reason, we conjecture the following regularity property: If 
s = Ll is a C 1 domain, A £ C(Ll) and z £ L°°(Ll), then 


(3.8) 


u g C°' e (Ll) 
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for every 9 < 1. We now derive the following regularity properties for the fractional 
optimal state u and the optimal adjoint state p. To accomplish this task, we define 


(3.9) A(u, p, Ud) — ||u|| H <.(n) + ||p||iP(n) + || z IU“(n) + IILid||z,°°(fi)- 

Theorem 3.6 (boundary regularity of u and p). Let u d £ L°°(Ll). If s £ (0, |), fl 
is a C 1 domain and A £ C(f2), then both u and p belong to C'°’ 2s (f2). In addition, 

(3-10) [u]c°> 2 »(n) + [p]c°. 2 ' , (p) < A(u, P> u d)- 

On the other hand, if s £ (|, 1), LI is a C 1,2 **- 1 domain and A £ C'°’ 2s_1 (ri), then 
both u and p belong to C 1,2s_1 (f2). In addition, 


(3-11) [u]ct. 2—i(fi) + [pIc 1 - 25 - 1 ^) < ^(u, p, Ud). 

In both inequalities the hidden constant is independent of the optimal variables. 


Proof. Let s £ (0, i). The results of Lemma 3.5 imply that u £ C' 0,2s (fi) and that 


(3.7) hold s. In view of the fact that u — Ud € L°°( fi), we apply, again, the results 
of Lemma 3.5 and conclude that p £ C' 0,2s (fl), and 

[p]c 0 ' 2s (O) < l|P||H«(fi) + INU“(n) + l|Ud||L“(n)- 


This estimate, in view of (|3.7[ ), allows us to derive (3.10). Analogous arguments 
can be applied to obtain (3.1l| for s £ (|, 1). □ 


We now present the following improved regularity result for the optimal control z, 
which shows that, for s £ (|, 1), and under some suitable assumptions, z £ C 0,1 ^). 

Theorem 3.7 (boundary regularity of z). Let u d £ L°°(Ll). If s £ (0, |), LI is a 
C 1 domain and A £ C(Ll), then z £ C' 0,2s (fi) and 


[ z ]c , °> 2 »(f2) < ||p||H»(n) + IH|ioo( n) + ||u d ||L=o(n). 

On the other hand, if s £ (|, 1), LI is a C 1,2s_1 domain and A £ C' 0,2s-1 (fi), then 
z £ C 0,1 (f2) and 


■Jc°> 


L (h) £ 


P||h»(0) + H u llL“(f2) + ll u d||L“(n)- 

Proof. The desired regularity results follow from Theorem |3.6| in conjunction with 

(fi) into C°’ a (Ll) 

□ 


the fact that the projection formula (3.6) maps continuously C 0, ‘ 
for er e (0,1]. 


We conclude this section with a conjecture regarding the regularity of the optimal 
control z when s = |. If £ L°°(Lt), Li is a C 1 domain and A £ C(f2), then 

(3.12) z £ C°’ e (Ll) 

for every 9 < 1. 


4. The truncated optimal control problem 


The state equation (3.2) of the extended optimal control problem is posed on 
the semi-infinite domain C = LI x (0,oo). Therefore, it cannot be directly approx¬ 
imated with finite-element-like techniques. However, since the solution to (3.2) 
decays exponentially in y Si Proposition 3.1], by truncating C to Cy = Ll x (0, 0 r ), 
for a suitable truncation parameter <y > 1, and setting a homogeneous Dirichlet 
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condition on y = y, we only incur in an exponentially small error in terms of y [5] 
Lemma 4.6]. We briefly review the results of [5] Section 4]. To do this, we define 

H 1 L {y a ,Cy) ■= {u £ H 1 (y a ,Cy) : v = 0 on d L Cy U 12 x {y}} . 

Then, the truncated optimal control problem reads as follows: Find min J(ti'QV, r) 
subject to the truncated state equation 

(4.1) ay{v,4>) = (r,trn0) H -, (n) xH«(fi) V(/> £ & l L {jj a ,Cy) 

and the control constraints r £ Z ac i. The bilinear form ay is defined by 

(4.2) ay{w,<j>) = j- f y a A(x)Vw • V</> Vie, <f> £ H^(y a ,Cy). 

0‘S J Cy 


We define the truncated control-to-state operator H : H _s (12) 3 r i-> trn v G 
H s (12), where v = v(r) denotes the unique solution to (4.1|. The map H is linear 
and continuous; see m Proposition 2.1]. With this operator at hand, we define 
the reduced cost functional j ■ Z ad 3 r -» j( r) G K by 


(4.3) 


j(r) = J(r,Hr), 


which is continuous and convex. In addition, the quadratic structure of j implies 
that j"(q)(r, r) does not depend on q and is positive definite, that is 


/'(q)(r, r) > i?||r||| 2 (n) Vr £ L 2 (12). 

O 1 

The truncated adjoint state p = p( r) £ Hf(y a ,Cy), associated with the extended 
state v = u(r), is defined as the unique solution to 


(4.4) ay(<j>,p) = (trn v — u d ,tro </>)i 2 (n) V0 £ Hl(y a ,Cy). 

We then have the following result [5] Theorem 4.5]. 


Theorem 4.1 (existence, uniqueness and optimality system). The truncated opti¬ 
mal control problem has a unique solution (v, r). The optimality system 


(4.5) 


'v = v(r) £ Hy(y a ,Cy) solution to (4.1), 


p = p(r) £ H\(y a ,Cy) solution to ( |4.4| ), 
f£Z ad, (tr fi p + dr, r - r) L 2 (n) > 0 Vr £ Z ad , 


holds. These conditions are necessary and sufficient. 


The following exponential approximation properties follow from [5} Lemma 4.6]. 

Lemma 4.2 (exponential convergence). If (%, z) £ H^(y a ,C) x HP (12) and (v,r) £ 
Hffiy a ,C.y) x (12) solve the extended and truncated optimal control problems, 
respectively, then 

P ~ ? IU 2 (0) < g-vWi^/ 4 (||r|| L 2 (n ) + ||Ud|| Z ,2(Q)) , 

IIV [ffi/ - v) || L 2 (yc t c) < e - ^^ 4 (||r|| L 2 (n) + ||u d || i 2 (n) ), 
where Ai denotes the first eigenvalue of the operator C. 


We conclude this section with the following regularity result for the truncated 
optimal control r and the truncated optimal adjoint state tr^ p. 
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Proposition 4.1 (Sobolev-regularity of r and trnp). Let r £ Z a( j be the truncated 
optimal control, be a convex domain and A £ C 0,1 (f2). If a < 0 < b and 
u d £ H 1_s (n), then r £ Hq(Q). If, in addition, £ Hq(Q), then trnp £ H K (fi), 
where n = min{l + 2s, 2}. 

Proof. The techniques of [41] Remark 25] allow us to transfer the regularity results 


of Lemmas |3.3| and [3T4| to r and trnp, respectively. To elucidate these results, we 
write, in view of separation of variables, a representation formula for the solution 
to pTj): v(x',y) = Yjk v kVk{x')xk{y ), where \k solves 


(4.6) 


Xk 


+ ay 1 x'k - AfeXfe = o, Xfc(0) = 1, Xk{y) = o, 


and {ipk} denote the eigenfunctions of the operator C. If I s and K s denote the 
modified Bessel functions of first and second kind [2 Section 9.6], then 

Xk{y) = (ak, s K s (-\f\ky) + b k , s I s {\/^ky)^ ■ 

The arguments of SD Remark 25] thus reveal that bk, s = — c s K s (^Xk‘y)I s (\f\k‘y) 1 
and that a,k, s = c s = 2 1-S /T(s). Define ek, s = 2 1_s 6fe !S /T(s). Then, in view of the 
properties of the Bessel functions we conclude that {e/c jS } decays exponentially to 
0 as k t oo, and in addition, that 

OO 

-lim y a v v (x',y) = V' A s k (d s - e k ,s}vkp(x'), 

yi0 k^i 


where d s = 2“T(s)/r(l — s). Since v solves problem (4.1), these arguments show 
that trn v solves the following fractional PDE: 

(4.7) trn v = r, 

where, for w £ Cg°(fi) and s £ (0,1), the fractional operator is defined by 

OO 

£ s w = ^ \ s k Wk(l - efc^dj 1 ); 
k= 1 


compare with (2.2). Since {efc jS } decays exponentially to 0 as k f oo, in particular 
is uniformly bounded, and then, we can extend the operator £ s to the space H S (D) 
defined by p3| . 

Lemma 3.5 in 15 


We can thus proceed, on the basis of (|4.7|) , as in the proof of 
to derive that 


r £ iLg(Sl). In view of this result, the arguments 
of Lemma 3.4 can be applied to derive that trnp £ H K (f2) for k = min{l + 2s,s}. 
For brevity we skip the details. □ 

The regularity properties derived by Caffarelli and Stinga in m for u, the solu¬ 
tion to (1.3), are based on the the fact that u = trn ft? and that ^ solves (3.21; see 
the proof of in Theorem 6.1]. In fact, since these arguments rely on the sctructure 
of the operator involved in the extension problem and are local in the extended di¬ 
mension, they can also be applied to derive regularity properties for trn v. where 
v solves (4.1). This, on the basis of Theorem 3.7 provides the following Holder 
regularity results for the truncated optimal control r. 

Proposition 4.2 (Holder-regularity of r). Let r £ Z a( j be the truncated optimal 
control. Let u d £ L°°(H), H be a C 1 domain and A £ C(H). In addition, if for 


s £ (|, 1), we have that H is a C 1,2s 1 domain and A £ C 0,2s 1 (H), then 


(4.8) 


r £ <™(H) for s £ 0, 


1 


and r £ C 0,1 (f2) for s £ 


1 
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For s = |, we conjecture the following regularity result. If £ L°°(f2), flis a 
C 1 domain and A £ then 

(4.9) f£C'°’ e (fi) 

for every 6 < 1 . 


5. A FINITE ELEMENT METHOD FOR THE STATE EQUATION 


In the next section we will propose a fully discrete scheme to approximate the 
solution to the optimal control problem (1.2 1 (1.4). The analysis relies, first, on 


the localization results of Section [3j and second, on finite element approximation 
techniques for solving (4.1) on curved domains; the latter being an extension of 
the priori error analysis developed in ESI ■ We comment that such an analysis is 
not trivial, since involves anisotropic meshes in the extended dimension and the 
nonuniform weight y a (a = 1 — 2 s £ (— 1 , 1 )), which degenerates (s < 1 / 2 ) or blows 
up (s > 1 / 2 ). 

It is instructive to review the results of [3D], which assume that fl is a convex 
polytopal subset of R" (n > 1) with boundary dfl. To do this, we start by recalling 


the regularity properties of and v, solutions to (3.2) and (4.1), respectively. The 
second order regularity of is much worse in the extended direction. In fact [401 
Theorem 2.7] (see [3H Remark 25] for v) yields 


(5-1) ||A a; /^|| i 2(j / a, jC ) + KV X ^|| W , C ) ^ ||/||Hl-=(fi), 

(5-2) \\%y\\L^ fi) ^ ll/lli 2 (n), 

with /3 > 2a + 1. These regularity estimates have important consequences in 


< 


the design of efficient numerical techniques to solve (4.1); they suggest that graded 
meshes in the extended (n+ l)-dimension must be used. We recall the construction 
of the family of meshes {SXy} over Cy used in [51l3Tj. First, we consider a partition 
Xy of the interval [ 0 , y\ with mesh points 


(5.3) 


y k = k 1 M 1 y i k = 0, 


,M, 


where 7 > 3/(1 — a) = 3/(2 s) > 1. Second, we consider = {K} to be a 
conforming mesh of SI, where K C M n is an element that is isoparametrically 
equivalent either to the unit cube [0,1]" or the unit simplex in R". We denote 
by Tq the collection of all conforming refinements of an original mesh We 

assume that Tn is shape regular m- We then construct a mesh STy over Cy as the 
tensor product triangulation of SAi £ To and Xy. We denote by T the set of all the 
meshes obtained with this procedure, and recall that T satisfies the following weak 
shape regularity condition: If Tj = K\ x Ii and T 2 = K 2 x I 2 £ STy have nonempty 
intersection, then there exists a positive constant <jy such that 

( 5 . 4 ) h h hj^<ay, 

where h[ = |/|. This weak shape regularity condition allows for anisotropy in the 
extended variable y mmm- 

Given a mesh £ To, we denote by lA/(5n) the set of all its interior nodes. 
We define hyr n = ma xxeSh k-K- For STy £ T, we define the finite element space 

( 5 . 5 ) V{& r ) = {W £ C°(Cy ) : W\ T £ Pi{K) ® Pi(/) VT £ W|r D = 0} , 
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where U 12 x {y} is the Dirichlet boundary. The set V\(K) is Pi (if) - 

the space of polynomials of total degree at most 1 - when the base K of an element 
T = K x I is a simplex. If K is a cube, ’Pi(A') stand for Qi(AT) - the space of 
polynomials of degree not larger than 1 in each variable. We also define the finite 
element space U(i5fo) = tro Y(^). We assume that ~ M n . This, in view of 
the fact that ff££y = M#=%, implies that ff£Ly ~ M n+1 . 

The Galerkin approximation of (|4.1|) is the function V £ V(A fr y) that satisfies 


(5.6) 


ay(V, W) = (r, tr n W) L 2 (Q) VW € V(^), 


where ay is defined in (4.2). We present the a priori error estimates derived in 
Theorem 5.4] and 0U| Corollary 7.11]. 


Theorem 5.1 (a priori error estimates). Let Sfy £ T and V(t3y) be defined by 
(5.5). If^( r) £ H\(y a ,C) solves (5.6) with z replaced by r, then 

(5.7) \\V(W(r)-V)\\ L , iyafi) < |log(#^)| fl (#^ r )- 1 /( n+1 )||r|| H i-. (n) , 


where y ~ log {ffLTy). Alternatively, if u(r) denotes the solution to (1.3) with 
forcing term r, then 

||u(r)-tr n V|| H . (n) < |log(#^)| s (#^)- 1/(ri+1) ||r|| H1 -»(n). 


Remark 5.2 (domain and data regularity). The results of Theorem 5.1 hold only if 
r £ H 1_s (12) and the domain 12 is sufficiently regular, for instance, convex. 


5.1. A priori error analysis for fractional diffusion on curved domains. 


In order to guarantee the regularity results of Theorem 3.7 and Proposition 4.2 


we 


need the following smoothness assumptions on the domain 12 and the matrix A that 


defines the operator £ in (1.5): 


(a) For s £ (0, g), 12 is a convex C 1 domain and A £ C(fl). 

(b) For s £ (i, 1), 12 is a convex C b 2s_1 domain and A £ C' 0 ’ 2s_1 (12). 


For the critical case s = 1/2, we assume (4.9), which in turns requires that Q is 
a convex C 1 domain and A £ Since we will be working on the basis of 

assumptions ([a]) and (Jbj) we cannot consider the domain 12 to be a convex polytopal 
domain in R" (n> 1). Instead, we consider a family of open, bounded and convex 
polytopal domains {12^}, based on a family of shape regular triangulations {^h}, 
made of simplices, that approximate 12 in the following sense: 


(5.8) 


C 12^, n <912^ C <912, |12\12^ 




where denotes the set of all the nodes of the mesh we refer the reader 

to |T8] or 02 Chapter 5.2] for details; see also 09]. From now on, we assume 12 to 
be a convex C 2 domain; the convexity property implies that {12^} C 12. 

Remark 5.3 (previous results and regularity of 12). We remark that, since the error 
estimates derived in 0] are based on the 77 1 -regularity of the optimal control z, the 
previous construction of the sequence {12^} is not needed in 0]: 12 can be taken 
as a convex polytopal domain in K™ (see Lemma 3.3). In contrast, we will operate 


under the regularity results of Theorem |3.7| and Proposition |4.2| and therefore, in 
order to have the validity of such results and handle the curved domain, we assume 
12 to be a convex C 2 domain. 
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On the basis of the previous construction, it is thus necessary to modify the 
definition of the finite element space Y(f7y). For the sake of simplicity, we keep the 
notation and define 


(5.9) V(^ y ) = {IV G C°(C» : W\t € Pi (K) ® P,(7) VT, ^|o\ n ^x(o,r] = °) > 

where Cyr = Llsr x (0, fk), T € ^ and the mesh STy of Cy is constructed as the 
tensor product of and Iy\ the latter being defined in Section [5j The discrete 
state equation then reads: Find V £ V(^) such that 


(5.10) 
where 

(5.11) 


a^(V,W) = (r,tr Q W) L2(a) VWeV(^y), 


a*r(V, W) = -j- [ y a A(x)X7V ■ VW. 

“s JCv 


We now present an extension of the a priori error estimate (5.7) of uni Theorem 
5.4]. In contrast to (5.7) the derived estimate allows us to consider curved domains. 

Lemma 5.4 (energy-error estimate on curved domains). Let f/(r) £ Hj^(y a ,C) 
be the solution to ( |3.2| ) with z replaced by r, and let V £ Y(37y) be the solution to 
(5.10). If Cl is a convex C 2 domain, A £ C'°’ 1 (fi), r £ L°°(Cl), and J ~ | log(#5^)|, 
then 


(5.12) ||V(^(r) - V)\\ L , {ya>c) < | log(#^ r )| s (#^ r )- 1 /( n+1) , 

where the hidden constant is independent of ^{r), V, r and Sfy. 


Proof. We start with an application of the triangle inequality and the exponential 
estimate of j40j Theorem 3.5] to deduce that 


(5.13) ||V(^(r) - V)\\ W , C) < ||V(^(r) - v)\\ LHyafi) 

+ l|V(« - V)\\ L , {yaM < e-^/ 4 ||r|| H - s(n) + \\V(v - V)\\ L , [v * fir) , 

where v corresponds to the the solution to (4.1) and Ai denotes the first eigenvalue 
of the operator C. It thus suffices to control the second term on the right hand side 
of the previous expression. To accomplish this task, we write 

||V(u - V)\\ L2{yafir) = ||V(v - V)\\ LHyaM + ||V(w - V)\\ L 2 {ya>eAc ^ 

and estimate each term separately. We begin with ||V(u — V r (r))|| i 2( y c, jC , J .), where, 
we recall that Cyr = Clyr x (0,ff). We denote by W e the extension by zero of 
W £ Y(3fy) to C. Then, given IT £ V(d£y), we invoke problem ( |5.10| and problem 
(4.1) with 4> = W e to arrive at the following Galerkin orthogonality property: 

asr(v — V, IT) = 0 VIT G V(^). 

This immediately yields 


||V(u - V)\\ L , {yaM = w inf ||V(u - W)\\ LHyaM , 

WeV(Sfy) 

which, in view of the piecewise polynomial interpolation results of EDI HU, and the 
regularity results of [401 Theorem 2.7], implies the following quasi-optimal error 
estimate in terms of degrees of freedom: 

(5.14) IIV(f - V)\\ LHyaM < I log7VrT- 1 /("+ 1 )||r|| H1 - 5(n) ; 

we refer the reader to the proof of uni Theorem 5.4] for details and remark that, 
the results of [101 Theorem 2.7] are valid under the assumption r £ H 1-S (r2). 
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We now bound ||V(u - V(r))\\ L 2 ( y c ^c y \C^)'- Since V(r) = 0 on f! \ Slj x (0,7], 
an application of the exponential estimate SSI Theorem 3.5] implies that 


(5.15) ||V(u - V)\\ L . {yafiy \Cr) = 

< e-^/ 4 ||r|| H -. (n) + ||W(r)|| L ^ fiv \C^y 

To control the remainder term, we use pointwise estimates for the harmonic exten¬ 
sion ^(r) that relies on the fact that r £ L°°(Q). These estimates are described, 
for instance, in Sj: |V a ,' £ ^'(r)| < y “for (x',y) £ Cy (SI inequality (6.1)]) and 
\d y ^(r)\ < y~ a for (, x',y ) £ Cy (SI inequality (6.2)]). These estimates imply that 


(5.16) ||W(r)||| 2(y<J)C \ Cjr) = ^ [ y a \V x ,W(r)\ 2 dx'dy 

Jo Jn\n°r 


o J r2\r2 op 

+ T [ tf a \d y W(r)\ 2 da/dy<7 1 - a \n\n J r\. 

J 0 J a? 

Since < h 2 $- and hyr n ~ N~ 1 ^ n+1 \ a collection of all the derived estimates 

allow us to conclude that 

l|V(^(r) - V)\\ L ^ yafi] < e-^/ 4 ||r|| H - s(n) + | log(#^)| s (#^)- 1 /("+ 1 ). 


This, in light of the fact that y ss | \og(ff7y)\, implies the desired estimate (5.12) 
and concludes the proof. □ 


Remark 5.5 (regularity of r). Examining the proof of Lemma 5.4 we realize that the 
critical step where the L°°(fl)-regularity of r is needed is (5.16). This assumption 


guarantees the pointwise estimates for ^(r) used to control its energy on Cy\Cyr- 


We present an improvement on Lemma |5.4| an error estimate that is quasi 
optimal in terms of approximation and only requires the H 1_,s (f2)-regularity of r. 
This improvement will allow us to derive an A 2 (n)-error estimate via a duality 
argument. 



Figure 1. An element A', for n = 2, such that K n dClyr ^ 0 and 
the corresponding curved region Ok- Sk denotes the side of AT 
that lies on dflyr, £k the corresponding arc formed by the curved 
boundary dfl and Ok the region bounded by E^- and Sk- 

Theorem 5.6 (energy-error estimate on curved domains). If r £ H 1-S (fl), then, 
under the framework of Lemma \ 5. 4\ we have that 

(5.17) ||V(^(r)-l/)|| i2( ^, c} <|log(#^)| s (#^)- 1 /("+ 1 )||r|| H1 - s(n) , 

where the hidden constant is independent of %(r), V, r, and 7y. 
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Proof. In view of the estimates (5.13)— ( [5TT5 ), we conclude that it suffices to bound 
the L 2 -weighted norm of V'^’(r) on Cy\Cyr- To accomplish this task, we proceed 
by a density argument that is inspired in the techniques developed in the proof of 
Lemma 5.2.3 in [33]. We present a proof for n = 2. Let K £ VLg- be a cell such that 
K n dflyr ^ 0. We denote by Sk the side of K that lies on the boundary dfl & and 
define Ok to be the domain bounded by the side Sk and the arc E k] see Figure]]] 
We choose local coordinates ((, y) such that £ is defined along the side Sk and y is 
orthogonal to Sk- We denote by rj = tp(£) the equation of the arc Sk- With this 
notation at hand, and for a smooth function <j>, we thus write the following relation 


<£(C,/b2/) = 0(C:¥>(O)J/) + / 

■M 0 

where y £ /, with I denoting an interval of the partition ly defined in Section [5j 
and (£, y) € Ok- Consequently, integrating on I and applying the Cauchy-Schwarz 
inequality, we arrive at 

[ y a 4‘ 2 ((,V,y)dy < [ y a <t> 2 {(,<p(0,y)dy + h 2 ^ [ [ y a (f) 2 (C,t,y)dtdy. 

Ji Ji JiJip(C) 

To obtain the previous estimate we have used the C 2 -regularity of the domain f2 to 
conclude that dist(x, dfl) for all x £ di}yr and then that |<£>(£) — y\ $h% 

Integrating, now, on Ok , we obtain that 


(5.18) 


ll < ^lli 2 ( 3 /°,OA'X/) ~ ^ 


m 2 L 2 (y a ,(dO K ndn)xI) 




x ,( P\\L 2 (y°-, 0 K xI)’ 


where we have used, again, that |</?(C) — y\ < /A,- Then, on the basis of a density 
argument, we obtain, for i = 1, • • • , n + 1, that 


(5.19) R,^(r)|| 2 < h%-Jd Xi nr)\\ L 2 (y a ,(50aHc?0) X/) 

+ L 2 (y a ,0 K xI)- 

We control ||9 a ;i‘2 f '(r)IU 2 (i/“ ) 9n X /) for i = 1,..., n. To accomplish this task, we 
use that H 1 (y a ,C.y) can be equivalently defined as the set of measurable functions 
w such that w £ H 1 {fl x (s,t)) for all 0 < s < t < y and for which the seminorm 
\\Vw\\ l 2 ( y a ,0?) is finite. Thus, if w £ H 1 (y a ,Cy ), an application of a standard trace 
inequality over a dyadic partition that covers the interval I implies 

[ V a w 2 <[ y a \\7w\ 2 . 

JdQxI JQxI 


Summing up (5.19) over all the cells K such that K 0 912^ ^ 0 and over all the 
cells / £ ly, using the previous estimate and applying the regularity estimate (5.1) 
(see [4(71 Theorem 2.7]) we obtain, for i = 1,..., n, the estimate 


(5.20) \\d x ^(r)\\ 2 L2{y ^ CAC ^ < h%Jd x yW(r)\\ 2 L2{yaM 

+ h%- n \\'V x , d Xi fy'(r)\\ 2 L2 (y^^Csr) ~ ~ (#^r) n+1 

In the previous inequality we have used that hyr n ~ (#>%) _1//n and that #S?y ss 
■ M k, M n+1 , where M denotes the degrees of freedom of the mesh ly. 

The estimate of the term \\d Xi ^{r)\[ L 2 ^ ya ^ dOKr[d Q^ XI ^ for i = n + 1 in ( |5.19 ), 
follows from the formula (|3.4[) and the asymptotic estimates for the modified Bessel 
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function of the second kind I\ s . In fact, [40] equations (2.28)—(2.29)] imply that 
= -c s \f\k (\AfcJ/) #s-i (\/Afcy) , -^-ipk{y) ~ y~ a K /2 as yf 0 +. 


For brevity we leave the details to the reader. In view of this obtained estimate 
and (5.201, we conclude, for i = 1,... ,n + 1, that 


||5xi^(r)|||2( y o~ ~ (#^r) 


"■+ 1 r 


BP-'lfJ)) 


which, in view of the estimates (5.13)—(5.151, yields the desired result (5.171. 
We now provide an L 2 (fl)-error estimate for tro(^(r) — V). 


□ 


Theorem 5.7 (L 2 (fl)-error estimate on curved domains). Let £ H\(y a ,C) 
be the solution to (3.2) with z replaced by r and let V be the solution to (5.10). If 
It is a convex C 2 domain, A £ C' 0 , 1 (f2), r £ H 1 -S (fl), and J ~ | log(#d£y)|, then 

(5.21) || tr n (^(r) - V)\\ LHn) < | log(#d^)| 2s (#^-)||r|| H1 _ S(0) , 

where the hidden constant is independent off/(r), V, r and Sfy. 

Proof. We start with an application of the triangle inequality and the trace estimate 
(2.7) combined with the exponential estimate E£B Theorem 3.5] to derive that 

|| tr n (^(r) - V)||i 2 (n) < || tr n (^(r) - v)|| L 2 (n) + ||tr n (u- fOIU 2 (n) 

< e _v/XTr/ 4 ||r|| H -» (n) + || trn (v - V)\\ L 2 {n) . 

To control the remainder term, we define e = tro(u — V ) and £ = v — V, and 
denote by *P ,y n : L' 2 (fl) -A U(5n) the standard L 2 -projection. With this notation 
at hand, we write ||e|| L 2 (n) < ||e - <P,%e|| L 2 ( n ) + || s P -% e|| L 2 (0 ^ :) . The control of the 
first term follows from interpolation on Sobolev spaces [29] Lemma 2.4], standard 
interpolation results, the trace estimate (|2.7|) and (|5.12|. In fact, 


||e-q3 5h e|| i 2 (n) < /£Je|| H .(n) < h^\\W£\\ L ^ M < | logiV| s iV-^||r|| H1 - 0(n) . 

It thus suffices to control the term Wty3r n e\\L 2 (Q. 5r )- We argue by duality and define 

d£Hl(y a ,Cy): ay{(f,e) = (fp^e, tr n 0 ) L 2 (n) £ Hl(y a ,Cy). 

Then, by setting cf> = £ and using that tro £ = e together with the definition of 
*P, we obtain that 


(5.22) 


l|q3 % e|| 2 2(Q) = ay{£,9) < ||V£ h Hy ^c r )\\V(e - ©)|| L ^ M , 


where 0 denotes the finite element approximation of 6 on the space V(£7y). Ap¬ 
plying Theorem 5.6 together with the regularity estimates ( |5.2[ ) we obtain 

l|V(0- ©)|| WM < |log(#^)| s (#^)- 1 /("+ 1 )||fP5he|| HI - 5( n ) . 

Since the family {Uh} is quasi-uniform, an inverse inequality implies the estimate 
ll < P^h e llHi-»(n) < h ^\\^sr n e \\L^{n)- This, combined with ( |5.22| ) and d5.17[ ), yields 

Il^^e||i2 (n) < | log(#^)| 2s (#^-)-^ 1 +^/^+ 1 )||^^e|U2 (0) ||r|| H1 -, (f2) , 
which implies (5.21|) and concludes the proof. □ 
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We conclude with the following error estimates for the approximation of problem 


(1.3) that follow from an application of Theorems 5.6 and 5.7 To state these results, 


we define the following approximation of the solution to (1.3) with z replaced by r: 

(5.23) U = tr n V, 


where V denotes the solution to the discrete problem (5.10). The discrete space for 
U is defined as U sr n '■= tro Y(^). 

Theorem 5.8 (error estimates for fractional diffusion). Let u(r) £ H s (f2) be the 


solution to problem (3.2) with z replaced by r, and let U £ U(5h) be its numerical 


approximation defined by (5.23). If Lt is a convex C 2 domain, A £ C' 0,1 (n) and 
y ~ | log(#^)|, then we have the following error estimates: 

(5.24) ||u(r) — Z7 ||h« ( £ 2) < |log(#^)| s (#^)- 1 / ( ” +1) ||r|| HI - 5(Q ), 
and 

(5.25) ||u(r) - U\\ LHn) < | log(#^)| 2s (#^)-( 1+s )/(" +1 )||r|| Hl - s(n) , 
where the hidden constants are independent oftf/[r), V, r and Sfy. 

6. A FULLY DISCRETE SCHEME FOR THE OPTIMAL CONTROL PROBLEM 

The results of previous sections are important in two aspects. First, we were 
able to replace the original fractional optimal control problem by an equivalent one 
that involves a local operator and is posed on the semi-infinite cylinder C. Then, 
we considered the truncated optimal control problem, that is posed on the bounded 
domain Cy , while just incurring in an exponentially small error in the process; see 
Lemma |4.2| Since in Section [5] we have proposed and analyzed a finite element 
discretization to approximate the solution to the truncated state equation (4.1) on 


curved domains, it remains to propose an efficient solution technique to solve the 
truncated optimal control problem. This is the content of this section. 

We assume that is a convex C 2 domain, and introduce a new fully-discrete 
scheme that is based on the approximation of the optimal control by piecewise 
linear functions on quasi-uniform meshes. This is in contrast to [5|, where the 
optimal control is discretized with piecewise constants functions. To approxi¬ 
mate the optimal state, we use the first-degree tensor product finite elements on 
anisotropic meshes described in Section |5.1| To be precise, the scheme reads as 
follows: min J(trn V, Z), subject to the discrete state equation 

(6.1) a^(V,W) = (Z,tr n W) L 2 (n) VfbeV(^), 

and the discrete control constraints Z £ Z a d(^h), where 

ZadW = Zad n {z £ L°°(fV): z £ C(n*),z \ K e Pi (k)vk, = 0} , 

where ayr is defined in ( 5.10| ), K £ AAi, the discrete space V(^y) is defined in ( |5.9| ), 
and Pi (A) corresponds to the space of polynomials of total degree at most 1. For 
convenience, we will refer to the problem previously defined as the fully-discrete 
optimal control problem. 

We denote by ( V,Z ) £ Y(^y) x Z a d(<%) the optimal pair solving the fully- 
discrete optimal control problem. If we define 


(6.2) 


U := tr n V, 
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we obtain a fully discrete approximation ([/, Z) £ U(^?h) x Z ad (<^o) of the optimal 
pair (u, z) £ HF(fl) x Z ad solving the fractional optimal control problem (1.21 (1.4). 


We recall that the finite element space U (f?o) is defined as U(^o) = trn V(^). 

Remark 6.1 (locality). The main advantage of the fully-discrete optimal control 
problem is that involves the local problem (6.1) as state equation. In addition, it 


can handle multi-dimensions easily and efficiently - a highly desirable feature. 


We define the discrete operator H sr y : Z a d(^n) D(5n) such that, for Z £ 
2w(=%), it associates tro V £ U(c%), where V = V{Z ) £ V(fZy) solves (6.1). 
With this operator at hand, we define the discrete and reduced cost functional 


(6.3) jsr v : Z ad -> M, j^(r) = J(r, H^r), 

which is continuous and convex. In addition, the second derivative j'f (q)(r, r) does 
not depend on q and is positive definite, i.e., 

(6-4) ^(q)(nr) > ||r||| 2(n) Vr £ L 2 (fl). 

This property will be important in the error analysis provided in Section |6.1| 

We also define the discrete adjoint state P = P(Z) £ V(5y) as the solution to 

(6.5) a s (P,W) = (tinV-u d ,tvnW) L 2 {n) MW £ V(^). 

We present the following result concerning the existence and uniqueness of the 
optimal control together with the first order necessary and sufficient optimality 
conditions for the fully discrete optimal control problem. 


Theorem 6.2 (existence, uniqueness and optimality system). The fully discrete 
optimal control problem has a unique solution (V,Z). The optimality system 


( 6 . 6 ) 


V = V(Z) £ V(iZy) solution of ), 

P = P(Z) £ V(fZy) solution of ), 

Z £ Zad(t^h), (trn P + i)Z, Z — 2)^2^ > 0 



VZ£ Z ad (3T), 


holds. These conditions are necessary and sufficient. 


In the next subsection, we will provide an a priori error analysis for the fully 
discrete optimal control problem. The analysis relies on the regularity properties 
of the optimal control r derived in Propositions |4.1| and |4.2| 


6.1. A priori er ror analysis for the fully discrete scheme. In view of the 
results of Lemma 


4.2 


to control the difference z — 2 in the L 2 (fl)-norm, it suffices 
to derive an a priori error bound for the term ||r — Z ||L 2 (n)- To accomplish this 
task, we introduce the Lagrange interpolation operator I^ n as follows: Given a 
mesh £ Tq and a function if £ C(£ty), we define 

ftrni’W) = ^2 ^ x ' e 

zeyx.%) 


where {( f > z } is the canonical basis of the finite element space U(c5o) = tro V(fZy); 
see 23 for details. A simple application of the triangle inequality yields 

(6-7) ||r- Z|| L 2 ( o) < ||f- /^ n r|| L 2 ( o) + \\I^ n r - 2\\ L 2 { p,y 

To estimate the terms that appear in the right hand side of the previous ex¬ 
pression, we consider a suitable partition of the mesh that is based on an 
assumption about the structure of the active sets We divide in three 
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subsets, that contain the active cells, inactive cells and cells with kinks , and are 
defined, respectively, as follows: 

( 6 . 8 ) = {K £ : r\ K = a or r\ K = b}, 

(6.9) &£ = {K e & n :a<r\ K < b}, 

(6.10) 3$ = \ (<% U ^ n 2 ). 

We notice that is the set of all the cells of the mesh that contain the 
free boundary between the active and the inactive set. We assume the following 
condition on the mesh and the optimal control r: there exists a constant C > 0, 
which is independent of hg- a , such that 

(6.11) ]T \K\<Ch^. 

This assumption is valid, for instance, if the boundary of the sets {x' £ fl : r(x') = 
a} and {a :' £ fl : r(x') = b} contain a finite number of rectifiable curves [571155 ,. 

We now bound the first term on the right hand side of ( |6.7| ). To simplify the 
exposition of the material, we define 


( 6 - 12 ) 


where 9 < 1. 


a := 


( 1 

2 ’ 

s e (o, i), 

2s, 


0 , 

S - I 
* 2 ’ 

.1, 

s £ (s> 1)> 


Lemma 6.3 (interpolation estimate). Let fl be a convex C 2 domain and A £ 
C°’\fl). If £ L°°(fl) D Hq^H) and a < 0 < b. then we have 


(6.13) 


|r-/%r|| L 2 (n) < h%. 


+cr 


where a is defined in (6.12), I^- a denotes the Lagrange interpolation operator and 
the hidden constant is independent of the continuous and discrete optimal variables 
and the mesh ST<y. 


Proof. We begin by invoking Proposition 4.2 to conclude that r £ C(fl) and then 


that Iy n r is well-defined. We now split the square of the error ||r — Isr^W l 2 (Q.) into 
two contributions: 


ll r — -^.% r lll 2 (n) — ll r — ^^h r lll 2 (n\n sr ) + ll r ~ -^% r ll! 2 (n^)- 
The first term on the right hand side of the previous estimate is bounded by 
using the fact that I,% :i I vanishes on fl \ fl^r and that r £ Z a d C L°°(fl ): 

(6.14) \\f ~ Isr u r\\%(n\ nsr ) 1$ 1^ \ ^ srWYAll-(n) < h\ max{|a|, |b|} 2 . 

We now control the remainder term ||r — I$cJ\\L 2 (n#y A first estimate can be 
derived in view of the results of Proposition |4.1| In fact, a standard interpolation 
estimate implies that 

(6-15) ||r ^a^\\ 2 L 2 {Qsr) ~ Mffi(O)- 

We remark that this step requires that a < 0 < b, in order to guarantee the H 1 
regularity of the optimal control r; see Proposition |4.1[ However, such an estimate 
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can be improved, for s G 1), by invoking the regularity results, on Sobolev and 
Holder spaces, of Propositions |d.l| and |4.2[ re spectively. To accomplish this task, 
we follow and invoke the partition (|6.8[)-(|6.10| of the mesh to write 


l r - J % r lli=(n er) 


= Y, H n ? - J ^ 


f! r lli 2 (K) — I ' 


II + III. 


J =1 


We proceed to control each term separately. We start with the term I: For each 
cell K G we have that r — = 0, consequently 1 = 0. 

If K G , we have that a < r\ K < b. Therefore r\ K = —{!/•&) tr^plx and then, 
over such a cell K , r and trap possess the same regularity, which is dictated by 
Proposition 4.1 ti 'np £ H K (f2), where k = min{l + 2s, 2}. We then arrive at 

trnp||i 


2At II 


n< E "k 




qn)- 


We now proceed to control t he te rm III. In view of the assumption (6.111, the 

4.2 and (4.9) imply, for s G (|, 1), that 


regularity results of Proposition 


hi < E i*iik-wi 

Kzsr* 


L°°(K) 


< h 


zr ll r llc°.‘ T (n) 


E 1*1 <h% 2 ' 


" r IIC7 0 -°-(O) ’ 


KeZT 3 


where a is defined in (6.12). We then collect the derived estimates for the terms I, 
II and III and invoke (6.15) to obtain that 

Ik ~ ~ 


upon realizing that k > ^ + cr for s G ( j, 1). This, in view of the estimate (6.14), 
implies the desired estimate. □ 

We now derive an instrumental result that will be important to estimate the 
second term on the right hand side of (6.7). To accomplish this task, we introduce 
the following auxiliary problem: 

(6.16) Q G V(5y) : ayr{Q, W) = (tr 0 u - u d , tr n IF) i 2 (n) MW G V(^), 


where v denotes the solution to problem ( |4.1| ) and ’V(lTy) is defined as in (5.9). 

Lemma 6.4 (auxiliary estimate I). Let Li be a convex C 2 domain and A G C 0,1 (f2). 
If Ud G H 1_S (H) and a < 0 < b, then we have the estimate 


(6.17) 


trn (p - P(r)) || ia ( n ) < | log(#^)| 2s (#^)-( 1+s V(«+D, 


where p = p( r) denotes the solution to (4.4) and the hidden constant is independent 
of the continuous and discrete optimal variables and the mesh SPy. 

Proof. We begin with a simple application of the triangle inequality to deduce that 

(6.18) || trf}(p- P(r))|| i 2 (n) < || tr n (p - Q)|| L 2 (n) + || tr n (Q - P(r))|| L 2 (n) . 

To bound the first term on the right hand side of the previous expression, we 
use the trace estimate (|2.7|) and the results of Theorem |5.7| We thus arrive at 


||tr n (p-Q)|U2( n) < ||V(p-Q)|| L 2 (ir , Cy) < | log(#^)H#^) 

We remark that, in view of the results of Lemma 
uniformly bounded. 


3.4 


• || — || 

" +1 1| trn u|| H i-»(n)- 

the term || tr n v||Hi-«(fi) is 
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We now proceed to bound the second term on the right hand side of (6.18). To 
accomplish this task, we invoke the trace estimate (2.7) and the stability of the 
problem that Q — P( r) solve to arrive at 

||trn(Q-P(r))|| L 2 (n) < ||V(Q - P(r))|| i 2 (j/ „ iCy) < || tr n (v - V(r))|| i 2 (fJ) . 

The remainder term is bounded by an application of of Theorem |5.7| This yields 


|| tr n (Q - P(r))|| L 2 (n) < | log(#^ r )| 2s (#^) »+i ||r|| H i-. (n) ; 

the H 1 _s (f 2 )-norm of r being uniformly bounded is a consequence of the fact that 
r £ Hq(SY); see proposition 


4.1 


Collecting the derived estimates, we arrive at (6.17) and conclude the proof. □ 

With the Lemm a |6.4| at hand, we proceed to estimate the second term on the 
right hand side of ( |6.7[ ). 

Lemma 6.5 (auxiliary estimate II). Let LI be a convex C 2 domain and A £ C 0 , 1 (f2). 
If Ud £ .ffo(^) Cl L°°(f2) and a < 0 < b, then we have 

l 


(6.19) 


ll^r-Z|| i 2 (n) <|log(#^)r/ l 2 , 


+ C 


where a is defined in (6.12) and the hidden constant depends only on the problem 
data. 


Proof. We begin with an application of the coercivity property (6.4) of jar. For an 
arbitrary q £ L 2 (f 2 ), we have that 

ni*j- z\\h(n) < i'kMV^r-zjsni-z) 

( 6 - 20 ) = jV(J 5 hr)(*Shf -Z)- j'-(Z){Ir a r - Z). 


We now invoke the variational inequality of the optimality system (4.5) with 
r = Z £ Z ad to derive that 

0 < -j'(r)(r-Z) = -j'(r)(Igr a r- Z) -/(r)(r- /%r). 

On the other hand, by setting z = Isr a r £ ^ad(=%) in the variational inequality 
of the discrete optimality system (| 6 . 6 |), we arrive at 


-j'^(Z)(I^r-Z)<0. 

In light of the previous two estimates, we proceed to control the right hand side 
of (6.20). In fact, we have that 

ni* a r - Z\\l* m < J^(/%r)(/ % r - Z) - - Z) - j'(r)(r - /%?) 

< - Z) - j^(r)(/ % r - Z)\ + [j' Sy { r)(i> n r - Z) 

- j'm^r - Z)\ - j'(r)(r - 1%r) =: I + II - III. 

Then, it suffices to estimate the terms I, II and III. To bound I, we invoke the 
Lipschitz continuity property of , which follows from the quadratic structure of 
. This yields, for p £ L 2 (fl), the error estimate 

|i| = li^W(W-bW-^)l 

(6.21) = |(H 5 y(/%r - r),H^(/_%?- Z)) L 2 (n ) + ^{Ig- n r - r, I- Z)l 2 (n)| 

^ \\IftJ~ r||L 2 (f 2 )||/%i ; - ^Hi 2 (n) < h a £ u 2 \\Iy n r - Z\\ L 2 ^, 










A FIRST-DEGREE FEM FOR A CONTROL PROBLEM OF FRACTIONAL OPERATORS 23 


where we have used the L 2 (fl)-continuity of the discrete control-to-state operator 
and the error estimate (6.13|. 

We now control the term II. Invoking the definitions (4.31 and (6.3) of j and 
j respectively, we write 

|H| = |(tr n (P(r) -p),Isr n f- Z) L i (Qj) \ < || tr n (P(r) - p)|| L 2 (n) ||/ 5h r - Z|| L 2 (f1) , 
and then, by applying the result of Lemma |6.4[ we obtain that 
(6-22) |II| < llog)#^)! 2 ^!/^- Z\\ L2(n) . 

We hnally bound the term III. To accomplish this task, we write 
III = /(r) (?-/_%?) = (tr n jp + i?r,r-/ t% r) i 2 (fl) 

= (tr n p + i?r,r-/ t% r) i 2 ( n\n a . ) + E III(A). 

Kesr u 

Since J^- n r vanishes on If \ and |fi \ < h 2 y , we conclude that 

(6.23) |(trnp + 'dr,r) i 2 (nvv) | < h 2 ^{\\ tr n p|U°= ( n) + IMU~(fi))||r||z,~(n); 


4.2 


the uniform control of || tropH^m) and ||r||L=o(n) follows from Proposition 
We now estimate each term III(A') depending on the location of the cell I\ £ 
with respect to the partition defined by ( 6 . 8 ) ( 6 . 10 ). 

• A £ <5^ : In this situation, the cell K is active. Consequently, (r — /%r)|# = 0 
and then III (A') = 0. 

• As : Since the cell A is inactive , the optimality condition of the system 
(4.5) immediately yields (trap + i?r)|x = 0. Consequently, III(A) = 0. 

• As : A first estimate can be derived by invoking the A 2 -standard projection 
operator tp & a : A 2 (fl^) —> U(^5h) as follows: 

III(A) = (tr 0 p + df-^Jtr n p + i?r),r-/^r) L 2 (if) 

+ WbhfafiP + •dr),fp j7n r - I& n f) L *(K)- 

We thus invoke standard interpolation and inverse estimates, to conclude, for 
s £ ( 0 , 1 ), that 

IH(P') ^ h 2 sr n (l|trnP + ^r|| ff i (if) ||r|| H i (if) ) . 

We notice that this estimate would only provide linear orden of convergence for 
the left hand side of (6.19). In what follows we present an improvement of this 
error estimate for s £ (j, 1). We begin with the estimate 

|IH(A)| < \K\\\tv u p + m\\ L ^ {K) \\r - IsrJ\\ L ^{K)- 

We bound || trap + To do this, we invoke the fact that K £ ^, and 

then that there exists a point x' 0 £ K such that (trop+dr)^) = 0. This, in view 
of Proposition 4.2 the assumption (4.9) and the definition of the Hdlder-norm, 
implies that 


(6.24) 


trnp + i?r|| L oo (x) = ||tr n p + dr- (tr n p + dr)(ai , 0 )|| L <x J(if) 


< A%||trnp- 


On the other hand, applying Theorem |3. 7| and Proposition 4.2 again, we obtain 


Ir- I 


■% r l|,L“(if) ^5 hZr n \\r\\ C 0 ,v(K) 
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Consequently, in view of ( |6.23| ) and the derived estimated for III(A'), we derive 
(6.25) |III| <h%,+ |HI(iC)| <h%- n +h% Y, \ R \^ h X +1 - 

KeST 3 Kesr 3 


Finally, inserting the estimates (6.21), (6.22) and (6.251 in (6.20), we derive the 
desired inequality (6.19) upon realizing that 1 + s > ^ + a for all s £ (|, 1). This 
concludes the proof. □ 


Theorem 6.6 (fractional control problem: error estimate). Let z and Z be the 
optimal controls for the fractional optimal control problem — ( |1.4[ ) and the fully 
discrete optimal control problem, respectively. If LI is a convex C 2 domain, A £ 
C 0 , 1 (fi) and \i d £ L°°(Ll) fl H 1 _s (fl), then we have that 

P - Z\\ LHn) < | log(#^ r )| 2s (#^)^+P, 


where the hidden constant is independent of the continuous and discrete optimal 
variables and the mesh tTy. 


Proof. We combine the exponential convergence result of Lemma 4.2 with the es¬ 
timates (6.13) and ( |6.19[ ) to arrive at 

P - Z\ |i 2 (fi) < P - 7 IU 2 (fi) + p - Z\\ L 2 (Q) 


< 


+ \l 0 g(#&r)\ 2a h 


2+ a 
sr ■ 


A natural choice of y comes from equilibrating the two terms on the right hand 
side of the previous expression: J ss |log(#Py)|. This gives the desired estimate 
and concludes the proof. □ 


Remark 6.7 (error estimate for s £ (0,1)). For s £ (|, 1), the estimate of Theorem 
16.61 reads 

P - Z\\ LHn) < | log(#^)| 2s (#^.)-™. 

Since the family is quasi-uniform, ~ M n and ~ #=% -M, we have 
that h*r n ss ~ (#5y) _1 /( rl+1 ^. This implies that the previous estimate 

can be rewritten as 

p-^IMa)<|log(#^)| 2s 4n- 

which, up to the logarithmic term, corresponds to the well known error estimate 
for first-degree approximation of the optimal control GE123- In the particular case 
that s = 1/2, Theorem | 6 . 6 | provides a slightly deteriorated error estimate: 

Wz-ZWvRtt) < |log(#^)| 2 p#^)F^(5 +e ), 


for any 9 < 1. For s £ (4,4), the result of Theorem 


6.6 


reads 


|z - Z\\ L 2(Q.) < |log(#^)| 28 (#^)^' 


> + 2s 


while for s £ (0, ^), Theorem 6.6 yields a linear order of convergence. 


Remark 6.8 (error estimates: n > 1). The advantage of deriving regularity proper¬ 
ties of the optimal controls z and r in Holder spaces has as a consequence that the 
provided a priori error analysis holds in a general n-dimensional setting, under the 
assumption that LI is a convex C 2 domain and A £ C , 0 ’ 1 (fi). 
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7. Numerical experiments 


In this section, we conduct a series of numerical experiments that illustrate the 
performance of the fully discrete scheme proposed in Section [ 6 ] and support our 
theoretical findings. 

The implementation has been carried out within the MATLAB software library 
«FEM j 15 . The stiffness matrices of the discrete systems (6.1 1 and (6.5) are as¬ 


sembled exactly, and the respective forcing boundary terms are computed by a 
quadrature formula which is exact for polynomials of degree 4. The resulting lin¬ 
ear system is solved by using the built-in direct solver of MATLAB. To solve the 
minimization problem, we use the projected Broyden-Fletcher-Goldfrab-Shanno 
(BFGS) method. The optimization algorithm is terminated when the ^ 2 -norm of 
the projected gradient is less than or equal to 10 -8 . 

To illustrate the error estimates of Theorem | 6 . 6 | we consider the following exact 
solution to the fractional optimal control problem (1.2) (1.4). Let n = 2, i? = 1, 
= (0, l) 2 , and A(x') = 1 in (1.5). The eigenvalues and eigenfunctions of the 
operator C are given by: 

= TT 2 (k 2 + l 2 ), y>k,l( x i, X 2 ) = sin(fc7nri) sin(Z7ra;2), 1,1 SN. 

To construct an exact solution, we consider the following modification or problem 
(1.3). Given s £ (0,1), the forcing term f and the control z, the aforementioned 
modified problem reads: Find u such that C s u = f + z in Q and u = 0 on dCl. If 
a = —0.5 and b = 0.5, and u,i = (1 + i?A 2 2 ) sin(27nri) sin(27r:E2), then we have that 
u = sin(27r.Ti) sin(27nE 2 ), p = — i?sin(27rxi) sin(27ra;2), and 

z = min {0.5, max {—0.5,— p/f?}} , 

where f = A 22 sin(27nri) sin(27rai 2 ) — z. We notice that z £ Hq(D.) D C' 0,1 (fi) for 
all the values of s £ (0,1). We remark that under this regularity property of the 


optimal control, the arguments developed in the proof of the estimates (6.13) and 


(6.19) guarantee the error estimate 

p-^||L 2 (0)<|l0g(#^)| 2s (#^ r )- 


2(i + l) . 


7.1. Piecewise linear versus piecewise constant approximation. In this sec¬ 
tion we explore the advantages of the fully discrete scheme of Section[ 6 ]when solving 
the fractional optimal control problem; we compare the performance of this scheme 
(Pi-scheme) with that of the numerical technique investigated in |5| that is based 
on piecewise constant approximation of the optimal control (P 0 -scheme). Table [l] 
show, for different meshes ^y, the error E z in the control approximation due to 
the Pi-schenre and the Po-scheme. Two values of the parameter s are considered: 
s = 0.2 and s = 0.8. #DOFs denotes the number of degrees of freedom of Sfy. It 
can be observed that, for a mesh Sfy with # DOFs = 137376, the error obtained 
with the Pi-scheme is almost an order in magnitude smaller than the corresponding 
error due to the Po-scheme of [5], 

7.2. Computational convergence rates for s £ (0,1). In Figure [2j we show 
the the asymptotic relation 

||z - Z\\ L 2 (n) K, (#&y)~ 5 

which illustrate the decay rate of our fully-discrete scheme of Section[ 6 ]for n = 2 and 
all the choices of the parameter s considered: s = 0.2, s = 0.4, s = 0.6, and s = 0.8. 
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# DOFs 

F Z (P 0 ;0.2) 

F z (Pi;0.2) 

F z (P o ;0.8) 

F z (Pi;0.8) 

432 

0.147712126 

0.131130828 

0.1482301425 

0.1470944750 

3146 

0.083305924 

0.036668665 

0.0840901319 

0.0443202090 

10496 

0.058953277 

0.020712242 

0.0588454408 

0.0241526956 

25137 

0.044253527 

0.012937511 

0.0441539905 

0.0148381456 

49348 

0.035650434 

0.008967500 

0.0356800357 

0.0101409325 

85529 

0.029769320 

0.007334747 

0.0297507072 

0.0080907623 

137376 

0.025419044 

0.005094037 

0.0254259814 

0.0056585074 


Table 1. Experimental errors for both: the fully discrete scheme 
studied in Section [6^1?! -scheme) and the P 0 -scheme proposed in 
[S]. Two different values of the parameter are s are considered: 0.2 
and 0.8. 


For s = 0.6 and s = 0.8, the presented results are in agreement with the error 
estimate of Theorem 6.6 For s = 0.2 and s = 0.4 the results of Figure [2] present an 
experimental convergence rate that is better than the one derived in Theorem |6.6| 
This is due to the fact, in this case, the exact optimal c ontr ol z £ Hg(f2) nC°’ 1 (fl); 
regularity property that is not provided by Theorem 3.7 for s £ (0, |). 


results are shown in Figure [3] for the following values of the parameter s: 
s = 0.3, s = 0.5, s = 0.7 and s = 0.9. 


Similar 
s = 0.1, 


\\z-Z\\v(a) 



Figure 2. Computational rates of convergence for the fully dis¬ 
crete scheme proposed in Section [6] on anisotropic meshes for n = 2 
and s = 0.2, s = 0.4, s = 0.6 and s = 0.8. The figure shows the 
decrease of the L 2 -norm of the error for the optimal control with 
respect to #d£y. In all the cases we recover the rate (#^) -1 ^ 2 . 
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||Z"-^IU a (£l) 



Degrees of Freedom (DOFs) 


Figure 3. Computational rates of convergence for the fully dis¬ 
crete scheme proposed in Section [6] on anisotropic meshes for n = 2 
and s = 0.1, s = 0.3, s = 0.5, s = 0.7 and s = 0.9. The figure 
shows the decrease of the L 2 -norm of the error for the optimal 
control with respect to In all the cases we recover the rate 

(#^)- 1/2 . 
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